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ABSTRACT 

Context. Accretion onto galactic and supermassive black holes occurs in different modes, which are documented in hard and soft 
spectral states, commonly attributed to an advection-dominated flow (ADAF) inside a truncated disk and standard disk accretion, 
respectively. At the times of spectral transition an intermediate state is observed, for which the accretion flow pattern is still unclear. 
Aims. We analyze the geometry of the accretion flow when the mass flow rate in the disk decreases (soft/hard transition) and evapo- 
ration of gas into the coronal flow leads to disk truncation. 

Methods. We evaluate the physics of an advection-dominated flow affected by thermal conduction to a cool accretion disk underneath. 
Results. We find re-condensation of gas from the ADAF into the underlying inner disk at distances from the black hole and at rates, 
that depend on the properties of the hot ADAF and vary with the mass accretion rate. This sustains an inner disk for longer than a 
viscous decay time after the spectral transition occurred, in accordance with the spectra that indicate cool gas in the neighborhood 
of the accreting black hole. The model allows us to understand why Cyg X-1 does not show hysteresis in the spectral state transition 
luminosity that is commonly observed for X-ray transient sources. 

Conclusions. Our results shed new light on the complex mass flow pattern during spectral state transition. 
Key words. Accretion, accretion disks - black hole physics - X-rays: binaries - stars: individual: Cyg X-1 



1. Introduction 

Accretion onto compact objects occurs in different modes. 
With the investigation of radiatively inefficient accretion flows 
(Narayan & Yi 1994), the observed spectra of low-mass X-ray 
binaries and Active Galactic Nuclei (AGN) could be understood, 
at low luminosities, as arising from a hot advection-dominated 
accretion flow (ADAF). The ADAF model or its variants pro- 
vide qualitative understanding of this hard spectrum from opti- 
cally thin emission (Narayan 2005). At high luminosities the ac- 
creting gas forms a standard accretion disk (Shakura & Sunyaev 
1973) reaching inward to the last stable orbit, producing a soft 
multi-color black body spectrum. 

But from the large number of observations for X-ray binaries 
during the last years (McClintock & Remillard 2006) it became 
obvious that for many sources the spectra document a highly 
complex appearance of accretion processes. Intermediate states 
were found, especially in connection with the change between 
the hard and soft spectral state (in both directions), or even the 
very high state, in both kind of sources, neutron star and black 
hole transients (Psaltis 2006). Hardness-intensity diagrams, e.g. 
flie diagram for the outburst of GX 339-4 in 2002/2003 (Belloni 
et al. 2005), show that the evolution of hardness with count rate 
is different during outburst rise and decay, indicating differences 
in the accretion flows via disk and ADAF. Similar changes in 
hardness are visible in diagrams for XTE J 1650-500 (Homan et 
al. 2003) and 4U 1630-47, (Tomsick et al. 2005). The variety 
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in diagrams for the hard/soft change can be seen from the re- 
cent compilation of Gierliiiski & Newton (2006). Kalemci et al. 
(2003) in their work on black hole X-ray transients during out- 
burst decay characterize the observed features as "several global 
patterns of evolution for spectral and temporal parameters be- 
fore, during, and after the transition". In general the Fe Ka line 
observed from stellar black holes and AGN, the reflection com- 
ponent, and also the timing properties might point to emission 
from cool gas in the innermost region or from the edge of the 
disk. 

Our work concerns the pattern of accretion flow at the time 
of spectral state transition. It is not yet clear how the accretion 
flow changes from one mode to the other, e.g. from disk accre- 
tion inward to the last stable orbit to a truncated disk with a hot 
advection-dominated flow in the inner region. We investigate this 
change starting from the picture that gas evaporated that from the 
disk feeds the coronal flow. What happens when the mass flow 
rate in the disk decreases (as during decline from an outburst) 
and becomes lower than the evaporation rate maximum? Then 
only a hot flow exists. We discuss that the disk reaching inward 
to the last stable orbit breaks up where the evaporation efficiency 
is maximal (at a distance of some hundred Schwarzschild radii 
{Rs = IGMjc^)). Farther inside, the cool disk would still ex- 
ist. however,But not being fed by a connection to the outside 
disk, such a relic inner disk would disappear in a viscous time. 
Apparently observed intermediate states last longer. An inner 
disk can then only survive if it is kept up by condensation of 
gas from the ADAF into the disk. Since these intermediate states 
are poorly understood at the present time we consider, as a first 
step, a new situation: Can cool gas exist underneath an ADAF? 
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We discuss the physics for this kind of accretion flow in an ana- 
lytical model. To include more detailed physics, an evaluation by 
numerical simulation has to be performed. We cannot yet answer 
in detail the question of how such cool gas could be recognized 
in the observed spectra. 

In the past the difficulty in interpreting the spectral obser- 
vations of transient low-mass X-ray binaries (LMXB) led to 
the question of whether the observations indicate "two indepen- 
dent simultaneous accretion flows", as discussed by Smith et al. 
(2002) and Pottschmidt et al. (2006). But whatever the geomet- 
rical configuration is in which accretion would happen simulta- 
neously via a hot flow and via a cool disk, at the same distance 
these flows would not be independent of each other. This es- 
pecially concerns thermal conduction. An interaction with cold 
clumps embedded in the hot flow causing thermal conduction 
of the ADAF in the radial direction was discussed by Yuan & 
Zdziarski (2005). 

One of the key questions in the study of the hot flow above 
cool gas is whether heat can be drained from the upper ADAF 
and be radiated away, so that matter condenses onto the cool 
disk. We investigate the physical properties of an ADAF affected 
by thermal conduction to an accretion disk below the ADAF. 
In particular, the energy balance of the hot flow at low height 
above the cool disk surface where ion and electron temperature 
couple is important. (For a discussion of the problem see Liu et 
al. 2006.) 

A peculiar aspect of such evaporation and condensation pro- 
cesses is the wide disparity between the amount of mass con- 
tained in the ADAF and in the disk, at each moment at each 
distance. This is a direct consequence of widely disparate tem- 
peratures in the two media. The velocities of the inward flow, 
proportional to the temperature, are thus also diff'erent by orders 
of magnitude and, for comparable mass flow rates lead to the 
very different column densities in ADAF and disk. How then can 
such a tenuous ADAF or corona, by condensation or evaporation 
affect and support a much more massive disk underneath? Even 
as a column of the ADAF loses only a fraction of its small mass 
during the short time that it passes over the disk surface, the un- 
derlying disk receives the contributions of many such columns 
before it has moved significantly inward, and the accumulated 
mass received can amount to a significant fraction of its own 
mass content. 

In Sect.2 we discuss the diversity of the outburst Ught curves 
of X-ray transients and the observations for an intermediate 
spectral state. In Sect. 3 we describe the accretion flow pattern 
related to the soft/hard spectral change. The physics of the hot 
flow affected by thermal conduction is treated in Sects.4 to 6. 
Using an analytical procedure we show in Sect.7 how the en- 
ergy balance in the radiating layer close to the disk surface de- 
termines whether evaporation or condensation occurs. In Sect. 
8 we determine the re-condensation rate of gas from the ADAF 
onto an inner reUc disk. The model allows us to evaluate how the 
condensation depends on the mass flow rate in the ADAF and 
the parameters involved. We draw conclusions regardingunder 
which circumstances condensation would lead to an inner weak 
disk below the ADAF, appearing in the spectrum as soft con- 
tribution. In Sect. 9 we show how our model can explain why 
Cyg X-1 does not have hysteresis in the transition luminosity, 
which is observed for all X-ray transient sources. In Sect. 10 we 
compare our model with observations for an intermediate state 
of X-ray transients, discuss the inner edge of the disk and the 
appUcation to disks in AGN. Conclusions are given in Sect. 1 1 . 

A better understanding of the pattern of the accretion flow 
should also give insight to the observed complex timing char- 



acteristics, the occurrence and type of quasi-periodic oscilla- 
tions (QPO) and the radio connections to the spectral states. Jet 
production corresponding to special accretion modes (GaUo & 
Fender 2005) also is part of the picture, arising from the accre- 
tion physics. 



2. The outburst light curves of X-ray transient 
sources 

2.1. The cause of the mass flow changes 

McClintock & Remillard (2006) examined 18 confirmed black 
hole binaries. In many sources the accretion rate seems to in- 
crease and decrease several times until finally after an outburst 
a quiescent phase follows. The large diversity of light curve fea- 
tures was documented by Chen et al. (1997), only a few sources 
showing the simple pattern of "fast rise - exponential decay" 
(FRED), but in many sources the luminosity varies in a complex 
way (see also the discussion by Gierliriski and Newton 2006). 
Two issues together create the variety of changes of luminos- 
ity and spectral states: (1) The outburst probably is caused by 
the ionization instabiUty, as the dwarf nova outbursts (Meyer & 
Meyer-Hofmeister 1984). The mass flow in the disk is related to 
the spread of heating and coohng waves initiating the changes 
between hot and cool disk structure. A larger disk size in X-ray 
transients together with irradiation leads to a more complex pat- 
tern in mass flow and outburst light curve (for the effect of a 
larger disk see accretion disks in symbiotic stars, Duschl 1986); 
(2) The mass flow rate in the disk resulting from this outburst be- 
havior is the key parameter that determines the accretion mode, 
for high rates via an accretion disk everywhere (soft spectrum), 
for lower rates a truncated outer disk and an ADAF inside (hard 
spectrum). 

2.2. The intermediate spectral state 

McClintock & Remillard (2006) redefined the X-ray states of 
black hole X-ray binaries in terms of quantitative criteria that uti- 
lize both X-ray energy spectra and power density spectra. They 
took as intermediate the spectra that did not belong to low/hard 
or high/soft state, nor to the very high state (sometimes a large 
part of the observation, e.g. 50% throughout the 2002-2004 out- 
burst of 4U 1630-47, Tomsick et al. 2005). We do not consider 
the intermediate state in the change to the very high/soft state. 
This latter state is of different physical nature (see the discussion 
by Fender et al. 2004, and the model of disk fragmentation by 
Meyer 2004). 

The observations show that generally the intermediate state 
phase during the hard/soft transition (in outburst rise) is shorter 
than that during the soft/hard transition (in outburst decay). For 
the former transition durations, of 3 to 11 days were found 
(Corbel et al. 2004, Gogu§ et al. 2004, Belloni et al. 2005). 

Our work focuses on the accretion flow pattern during the 
latter transition. Kalemci et al. (2003) analyze the evolution of 
spectral and temporal properties of several galactic black hole 
fi-ansients during outburst decay: XTE J1650-500, GRO J1655- 
40, XTE J1748-288, XTE J1755-324, 4U 1630-47, XTE J1550- 
564, XTE J1859-H226 and GX 339-4 (for this source see also 
Revnivtsev et al. 2000) and the spectral analysis shows a de- 
crease of the disk component until it becomes unobservable with 
PCA /RXTE within 15 days after the transition. Similar times 
were found for GRS 1758-258 (Smith et al. 2001, Pottschmidt 
et al. 2006). Cyg X-1 is a special case with an intermediate state 
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lasing for weeks (Cui et al. 1997, Zdziarski et al. 2002). The re- 
sults for neutron star sources show the same trend with shorter 
timescales (Barret & Olive 2002). 



3. The accretion flow pattern at the soft/hard 
spectral transition: a gap in the disk? 

In earlier work on accretion onto compact objects (Meyer & 
Meyer-Hofmeister 1994, Meyer et al. 2000, Liu et al. 2002, see 
also Rozanska & Czerny 2000a, 2000b) we found that above the 
cool accretion disk a hot corona always is built up, fed by mat- 
ter evaporating from the disk. The efficiency of evaporation, a 
siphon flow, depends on the distance from the compact star The 
evaporation rate increases with decreasing distance to the cen- 
ter until a maximum (Mevap)max is reached at /?evmax- The mass 
flow rate in the disk determines whether the coronal flow is im- 
portant. If it is higher than (Mevap)max the disk is essentially un- 
affected and reaches inward to the last stable orbit, and the hot 
flow is unimportant (this is the usual situation when the luminos- 
ity is high during an outburst, see also discussion in Sect. 8.3). If 
the mass flow rate in the disk is low, the disk is truncated where 
the evaporation has transfered all gas to the hot flow/ ADAF (the 
usual situation in quiescence). According to the mass flow rate, 
spectral transitions occur. Note that these two transitions occur 
at different luminosities, an observational feature known as hys- 
teresis in X-ray transient light curves. This is theoretically under- 
stood as due to the different hard/soft character of the irradiation 
from the central region in the two cases (Meyer-Hofmeister et 
al. 2005, Liu et al. 2005). 

When during the decrease from an outburst the mass flow 
rate in the disk becomes equal to the evaporation rate maximum, 
the disk truncation starts at the distance /^evmax, and a gap ap- 
pears, filled by an ADAF, as illustrated in Fig.[T]. With the further 
decrease of mass flow the gap should widen due to evaporation. 
What happens to the left-over inner cool disk? Diffusion causes 
a spread inwards and outwards, decreasing the amount of disk 
matter due to accretion onto the compact object and evaporation 
in outer regions. Only if matter can condense from the ADAF 
onto the disk can such a disk survive longer than the viscous 
time for these processes, which is only a few days in the inner 
disk. If the disks still exists the spectrum can show contributions 
from the ADAF and a reflection from the left-over inner cool 
disk. 



4. ADAF properties 

As our analysis of the hot ADAF above a left-over cool disk will 
show, the upper part, i.e. the main part in vertical extension, is 
not affected by thermal conduction to the cool disk below. We 
make use of the solutions of self-similar advection-dominated 
flows by Narayan & Yi (1995). The properties of these flows 
scale with black hole mass, mass flow rate and distance from 
the black hole. They further depend on the viscosity and the as- 
sumed magnetic field strength. Following Narayan et al. (1998) 
the magnetic pressure is written 



coronal flow 



(1) 



with /3 ratio of gas pressure to total pressure, p density and 
isothermal sound speed. Shearing box simulations of turbulence 
driven by the magneto-rotational instability in a collisionless 
plasma by Sharma et al. (2006) yield /3 values around 0.8. We 
use this value in our analysis. As the ratio of specific heats of the 
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Fig. 1. Accretion flow in soft, intermediate and hard spectral 
state (with decreasing mass flow rate). Upper panel: In the soft 
state accretion occurs via the disk, some gas is evaporated to a 
coronal flow, but re-condenses onto the disk; in the inner region 
only the disk flow exists, which causes the soft spectrum. Middle 
panel: In the intermediate state the ADAF and an inner disk con- 
tribute to the spectrum. Lower panel: In the hard state only the 
hot ADAF contributes to the spectrum. 



magnetized plasma we take y = (8 - 3y8)/(6 - 3yS) (Esin 1997) 
though the true value would require a more detailed analysis. 

For the chemical abundance a hydrogen mass fraction of 
0.75 was used. The solutions for pressure, electron number 
density, viscous dissipation of energy per unit volume and 
isothermal sound speed are (Narayan & Yi 1995) 



p = 1.87 X lO'V^'m^'mr^^'Vm" 



= 5.91 X W^a-'tn'mr 



(2) 



= 2.24 X 10^V"^mr"Vgscm"^s"\ 
cl = 1.67 X lO^^r^'cm^ s"^ 



where a is the viscous coefficient, m the black hole mass in units 
of solar mass M©, m the mass flow rate in units of the Eddington 
accretion rate Mndd = 1 -39 X lO'^m g/s, and r the radius in units 
of the Schwarzschild radius /?s - 2.95 x lO^m cm. The ion num- 
ber density is «, = n^lLOTIl . In an ADAF, ion and electron tem- 
peratures r, and Tf, closely follow 



1.98X lO'V'K 



(3) 



and, if is much smaller than T, this value can be taken for 
Ti alone. Besides viscous heating, compressive heating q'^ - 
jjT^q^ (Esin 1997) is important. 
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Fig. 2. The hot two-temperature ADAF together with the radi- 
ating one-temperature layer below the coupling interface. 

5. The electron thermal profile 

In a pure ADAF, at radii r < 10^, the coupling between ions 
and electrons becomes poor. Almost all of the accretion en- 
ergy stays with the ions, which reach near virial temperatures, 
the flow has a high vertical extent, and the density is low 
(see review by Narayan et al. 1998) . The electrons cooled 
by bremsstrahlung and synchrotron(-Compton) radiation stay at 
temperatures around lO'^-^K. Radiation losses are very weak. 

This situation changes if there is a disk below the ADAF. 
Due to the large temperature difference between the hot ADAF 
and the cool disk, thermal conduction results in further cooling 
of the electrons. While the ion temperature is not much affected, 
the electron temperature drops with height until, near the bottom, 
a temperature - T^pi is reached at which coupling between 
ions and electrons becomes efficient, and from then on ion and 
electron temperatures are the same, T, - - T, as illustrated 
in Fig.|2]. 

The electron temperature has practically no influence on the 
dynamics of the flow and the thermodynamics of the ion gas 
and is subrelativistic, kTf, < nifC^ (k Boltzmann constant, m,, 
electron mass, c speed of light). The rate of energy transfer from 
ions to electrons is given by Stepney (1983). Since in the two- 
temperature advection-dominated hot flow the ions are at a much 
higher temperature than the electrons a simplified formula (Liu 
et al. 2002) can be used 



qie = 3.59 X lO^^^^gcm^s^^deg"' n^n,T; 



= 1.05 X l0^^gcm-^s-^deg^'^T;^'^a-WWr-\ 



(4) 



to a first approximation the radiation losses and calculate the 
resulting thermal electron profile from the balance between ther- 
mal conduction and collisional heating alone. 



dF, 



dz 



= -qie(Te)- 



(6) 
(7) 



For the thermal conductivity coefficient kq we take the standard 
value (for the effect of /<•() on the coronal structure see Meyer 
& Meyer-Hofmeister 2006). Formally multiplying the left sides 
and the right sides of these two equations by each other and mul- 
tiplying by dz reduces this to a first order differential equation for 
FciJe) with the solution 



Fl = KQ{Kn,nJi){Tl 



"^^gcm^s~^deg^^^ 



(8) 



with K - 1.64 X 10" gem s~ deg ' . The integration constant 
T,„ is the maximal electron temperature reached with height. 
Note that the downward directed heat flow has a negative value. 
Integration of dzjdT^ (Eq|6l) using Eq.® yields the electron 
temperature profile in the ADAF, 



Zm Z ~ 



KriirigTi 



T 3/2 



•JTJT„, 



: dx. 



(9) 



For the low height z - Zcpi where has dropped to the cou- 
pling temperature Tcpi, in Eq.® we may neglect Zcpi compared 
to Zm, and with T^^JT^ small compared to one, the integral may 
be extended from zero to one to give 0.719. Eq.(|9]l then yields 
the maximal electron temperature in the ADAF, 



7^/2 = \39z„,{KninJilKQ) 



1/2 



(10) 



For the height Zm one may take the vertical scale height for 
which is 1/ V2 of the density scale height Cj/t^K (f^K Kepler 
angular velocity) of Narayan et al. (1998), 



1.27 X 10^ m rem. 



(11) 



With this value we obtain from Eq.lfTOl). Eq.®, neglecting 
r^pj compared to T^, finally yields the thermal flux drained from 
the ADAF that enters the coupling interface at z = Zcpi from 
above. 



F^"'^^ = -^f °(/:«,«,r,)^/'«(i.39zj 



2/5 



(12) 



Coupling of ions and electrons is reached when viscous and 
compressive heating of the ions is balanced by the collisional 
heat transfer from the ions to the electrons, qie - + (f . This 
then determines the coupling temperature r^pi. Since the cou- 
pling occurs at low height the density there is higher than the 
height averaged value (Eq. (|2|i), by a factor 2/ y[n (vertical den- 
sity distribution n - no exp(-z^ /H^)). This yields the coupling 
temperature as 



Fcpi = 1.98 X 10 V 



,9„-4/3^2/3 



= 1.24X WK, 



(5) 



the numerical value for q'=0.2 and m=0.02 (see Sect. 8.1). 

In the case considered here, conductive cooling limits the 
peak electron temperature in the ADAF to a value Tm some- 
what lower than the limit where radiation losses balance col- 
lisional heating. The synchrotron radiation is optically thick, 
which makes the synchrotron radiation losses drop steeply 
with temperature (Narayan & Yi 1995, Mahadevan 1997). As 
bremsstrahlung also decreases with temperature we may neglect 



6. The coupling interface 

In the case of no mass exchange between ADAF and disk, 
Eq.(fT2li gives the heat flux that enters at z = Zcpi into the ra- 
diating layer from the ADAF above. This is the borderline case 
considered by Liu et al. (2006). However if mass condenses from 
the ADAF onto the disk it carries the high thermal heat content 
of the ADAF ions to the interface. Energy conservation across 
the interface yields 



^ cpl 



r 



1 l%T 

■ r-i>\ y" 



7 



r ■ 



_i/M 

lySl J", I 



(13) 



The factor i arises when we include magnetic work and inter- 
nal energy advection together with the corresponding gas terms 
in the simplified ADAF parameterization of the magnetic field. 
Here Fcp\ is the heat flow leaving the interface on the lower side, 
is the gas constant, T is the common temperature of ions and 



F. Meyer, B.F. Liu and E. Meyer-Hofmeister: Re-condensation from an ADAF into an inner disk 



5 



log p 



logT 
11 



10 



Zcpi — ^ z 1 

Fig. 3. Changes of ion temperature, electron temperature and 
density at the distance of 80 Schwarzschild radii and for a mass 
flow rate of m=0.02 in the ADAF. The vertical height z is mea- 
sured in units of the vertical scale height for «,. 

electrons on that side, and m, is the mass flow rate of the con- 
densing gas per unit area. (It has the dimension of gcm^^s"' and 
should not be confused with the non-dimensional ADAF mass 
accretion rate m and similar quantities defined later) In conden- 
sation flow it has a negative value. The heat flow f cpi entering 
the radiating layer below the ADAF from above is thus increased 
over F^^^^, 



T7 17 

■fcpl - r, 



ADAF 



1 _ ^ 

m*. 



where 



y-1 
r 



^(1 - e) 



(14) 



(15) 



is a normalization mass flow rate per unit area (positive and of 
the same physical dimension as m,). It depends on the ADAF 
parameters. The ratio 



(16) 



is a small number, typically of the order of 10"^. An example for 
temperature and density distributions is shown in Fig. [3]. 

(In the case of evaporation (m. > 0) the situation at the 
boundary between the radiating layer and the ADAF is different. 
While for a descending flow, collisional coupling immediately 
heats the electrons and adds to their thermal flux, in an ascend- 
ing flow, on passing through the boundary, the ions becoming 
uncoupled must first be heated up by friction to the ADAF tem- 
perature before they are fully incorporated in the ADAF. Here 
we consider only the condensation case.) 

7. The radiating layer below the ADAF 

The parameters of the ADAF determine what happens in the 
layer between the coupling interface and the disk surface. 

Case (1): If the pressure is high, the density in the lower layer 
also is high. Then the conductive flux drained from the ADAF 
is already efficiently radiated away at some height before the 
disk surface is reached. Radiative cooling however continues and 
must be served by additional heat released from gas descending 
through the temperature profile to the disk surface, i.e. matter 
condenses from the ADAF onto the disk. 

Case (2): If the pressure is low, only part of the conductive 
flux can be radiated away; the remaining part is taken up by cool 



disk matter heating up as it rises through the temperature profile, 
i.e. matter evaporates from the disk into the corona. 

A borderline case in between, case (3), occurs if the ADAF 
pressure allows the thermal heat flux drained from the ADAF 
to be radiated away exactly on reaching the disk surface. Then 
there is neither condensation nor evaporation, i.e. no mass is ex- 
changed between disk and ADAF. 



7. 1. The energy balance in the radiating layer 

We derive the condition for the different cases from the energy 
balance in the radiating layer (compare Liu et al. 2002). We use 
a simplified form of the energy equation, keeping only the dom- 
inant contribution of internal heat, pressure work, and thermal 
conduction, together with bremsstrahlung cooling ne«,A(r), 



d_ 

dz 



m, h /< c 

'r- 1 1^ 



- -neHiKiT). 



(17) 



We express density by temperature and (constant) gas pressure 
P p. This value is taken at the bottom of the ADAF region, and, as 
for the density, is slightly higher than the vertical mean pressure 
(Eq.(|2]l). Assuming free-free radiation for > lO^^K, n(,n, A(r) 



becomes ^^{p pfbT-^'^ with b = 



j0-26.56g cm^s-^deg-'^^ 



(Sutherland & Dopita 1993). For simplicity we use this law also 
for smaller T . The justification is that because of the very steep 
temperature profile below lO^'^'K such regions contribute only 
a negligible amount to cooling of this layer (Liu et al. 1995). 
Contributions of gravitational energy release, frictional heating, 
and side- wise advection of mass and energy can be neglected for 
a small extent of this layer. Likewise, kinetic energy is negligible 
since at high density the flow is highly subsonic. 

To solve the second-order differential equation Eq.JTTIi. we 
use T as the independent variable and define a new dependent 
variable g{T) = KoT^^^dT /dz = -F^T . We now obtain the first- 
order differential equation 



2„2 



' d\nT 



0.25/3> 



■Kob ■ 
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'r - 1 



-ig - gi)(g - gi). 



with 



2 y- \ p 



2 y- \ p 



0250^ 



(18) 



(19) 



and g2 differing from gi only by the sign of the square root. 

This equation has to be solved with the upper and lower 
boundary conditions g - at T = Tcpi and F^ - Q al T - Q. 
In reality, the temperature does not drop to zero at the disk sur- 
face but the solution is practically independent of the exact val- 
ues at the lower boundary (except for a very narrow range in z at 
the bottom) as long as flux and temperature there become small 
compared to those at the upper boundary, with no consequence 
for our results, m. is the Eigenvalue to be determined. 

The only solution that fulfills the lower boundary condition 
is the singular solution g{T) = gi = const, i.e. a linear relation 
between -Fc and T. The ratio of these two quantities at the upper 
boundary determines the value of gi. This requires 



r-1 



_/rA 



1 -C- 



1 



1 -(2f)2 



(20) 
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with 



C - Kob 



[cpl_\ 



^2 ]\fADAFj 

Eqs.(l20]l. ( fTsT i, and ( fTSI l yield a quadratic equation for m 



(21) 



1- 



-e|l--i|-(l-e)C = 0, 



(22) 



with the solution 



m* 2 V 4 

The quantity C of Eq. (12111 compares the radiation loss in the ther- 
mal profile with the heat flux drained from the ADAF. For C = 1 
Eq.(|20]i gives m. - 0, the borderline case discussed by Liu et al. 
(2006). For C < 1, the radiation losses are too weak, m, is posi- 
tive and the heat flow is used up to heat the evaporating gas to the 
(un-)coupling boundary temperature r^pi. If, on the other hand, 
C > 1, the gas in this layer is efficiently cooled, sinks down and 
condenses. 



8. The mass exchange between ADAF and disk 

8. 1 . The borderline case: no mass exchange, C = 1 

The value of C depends on the ADAF parameters a, /3, m and 
r, but not on m. This shows that the same situation holds for su- 
permassive and stellar black holes alike. Also note a significant 
dependence on /3, C scales almost with (Eq. (l2T]) ). For y6=0.8 
we get 



C = 0.96a-28/i5^8/i5^-i/5. 



(24) 



The radius r - at which C - \ marks the transition between 
outside evaporation and inside condensation of an inner disk be- 
low the ADAF, 



= 0.815 a 



-28/3^8/3^ 



(25) 



Here, we take a=0.2 (corresponding to a—Q.'i in Shakura- 
Sunyaev notation) but note the extremely strong dependence of 
Yc on the choice of this parameter. Since a incorporates sev- 
eral underlying physical processes of unresolved detail its ex- 
act value is not known. Various values of the above order of 
magnitude are quoted in the literature (see Narayan et al.l998). 
The re-condensation process discussed here occurs in connec- 
tion with spectral state transitions for which Maccarone (2003) 
found m - 0.02 from a detailed investigation of these transitions 
in X-ray binaries. For these values C becomes 1 at a distance of 
about 80 Schwarzschild radii, an order of magnitude discussed 
for inner disks in intermediate states (Ibragimov et al. 2005). The 
dependence of C on mass flow rate and distance obtained allows 
us to derive conditions for the existence of an inner disk below 
the ADAF 



8.2. Condensation in the inner region. 

It is an interesting question as to how the fraction of the ADAF 
that re-condenses onto the inner disk and then flows inward via 
the disk depends on the ADAF parameters. For this we add all 
contributions from the distance Rc - rcRs where condensation 
begins down to an inner cut-off radius - r,7?s, e.g. the ra- 
dius of the innermost stable orbit around a non-rotating black 



(D 
T3 




0.05 



log r 



Fig. 4. Re-condensation of gas from the ADAF onto an inner 
disk as function of the distance r (in units of R5) for various mass 
flow rates m in the ADAF (in units of MEdd)- The curves display 
the spatial distribution of contributions of logarithmic intervals 
in r as areas below the curve. The ordinate is re-condensation 
rate relative to the mass flow rate in the ADAF. Re-condensation 
occurs at a larger distances for higher rii. Below 3Rs the curves 
are dotted. 



hole, r,- = 3. As discussed later, the inner edge may well lie 
at some larger distance. Like the ADAF accretion rate, the to- 
tal condensation rate will be measured in Eddington accretion 
rate. (In our definition condensation rates have a negative sign.) 
Using Eq.(|20li, (for e «c 1 ) the total rate of condensation onto 
both sides of the disk becomes 



^cond I ~ 



r''^ 4-nR^ dR 
I -: rhy — 

JR, MEdd ~ R 



-r 



.o(Vc-i)^ 



with 



47TR-m; 

MEdd 



= 8.31 X 10""q; 



3„-7/5™V/5 3/5 



m 



(26) 
(27) 

(28) 



The integrand yfc - 1 ) gives the contribution of logarithmic 
intervals in r. Its distribution is shown in Fig.|4]for various values 
of ih. The ordinate is the integrand relative to the ADAF mass 
flow rate ih. The area underneath each graph gives the relevant 
total condensation rate. 

Introducing x = r/r^ instead of r in Eq.(l26b and making use 
of Eqs. (l24l l. ( |25] ), and ( |28] l one can evaluate the integral and 
obtain the total re-condensation rate as 



m 



2.34 X IQ-^a-'^m^fiRJR^) 



with 

fix) = 1 - 6x^1^ + 5x^'^ . 



(29) 



(30) 



The function / approaches 1 when re-condensation begins at a 
radius R^ that is large compared to the inner cut-off" radius 
As the ADAF mass flow rate decreases R^ moves inward, / and 
the re-condensation rate, i.e. the mass flow rate in the inner disk 
decrease until at a certain lower ADAF mass flow rate the outer 
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disk edge Rc reaches Ri where / becomes and the disk size has 
shrunk to zero. For even lower mass flow rates there is no re- 
condensation and the inner disk disappears. This intriguing be- 
havior appears able to explain why Cyg X-1 does not show the 
hysteresis between hard/soft and soft/hard spectral state transi- 
tions that all transient X-ray sources display. We discuss this in 
the next section. 

8.3. Complete re-condensation? 

The model results depend on the parameters used to describe 
the ADAF. The dependence on the parameters a and m is par- 
ticularly strong for the re-condensation rate (Eq.(l29]l) and the 
radius at which re-condensation begins, i.e. the outer radius 
of an inner disk below the ADAF (Eq.dZSll). In the example 
given above both values are rather small. However, a small de- 
crease of a and a slight increase of ih can bring these values to 
100% re-condensation with re-condensation beginning already 
at 1000 Schwarzschild radii (this is an extrapolation as our inte- 
gral Eq.(l26ll neglects changes of m by re-condensation). 

This is relevant for standard disk evaporation models (e.g. 
Meyer et al. 2000). They predict that at distances of a few hun- 
dred to a thousand Schwarzschild radii there is a significant mass 
flow rate in an accretion disk corona, of a few percent of the 
Eddington value. If the mass accretion rate drops below that 
value, all the mass evaporates into the corona and from there on 
flows inward as a pure coronal flow, i.e. an ADAF. This yields a 
standard model for spectral state transitions. If the mass flow rate 
in the disk exceeds that critical value the disk is not truncated and 
continues to the last stable orbit. This accretion mode gives the 
characteristic soft state spectrum of a multi-temperature black 
body. But even in this case the rather strong coronal flow at the 
distances of a few hundred Schwarzschild radii must exist. Why 
does it not flow on as a 2-temperature ADAF above the disks and 
shows as a significant hard contribution to the soft state spec- 
tra? The answer from our intermediate state analysis is that re- 
condensation becomes complete for mass accretion rates as high 
as the coronal evaporation rate at a few hundred Schwarzschild 
radii. This explanation connects the processes of disk truncation 
and re-condensation. The upper part of Fig. 1 illustrates the sit- 
uation. 

8.4. Hysteresis 

Even for lower mass accretion rates, when the outer disk already 
is truncated, re-condensation onto the remaining inner disk may 
be close to complete and still yield an apparent soft state in the 
innermost region. This extends the soft state to a lower mass 
accretion rate on the declining part of an outburst and leads to 
lower luminosity for the transition from the soft to the hard spec- 
trum in decline than for the opposite transition during the rise, 
i.e. hysteresis. Such hysteresis would be independent of, and add 
to the irradiation caused hysteresis referred to in Sect. 9. 

9. Cyg X-1, a black hole accretor without hysteresis 

Cyg X-1 shows only moderate changes in the mass accretion 
rate, attributed to the fact that it is a wind accretor. This is very 
different from transient X-ray sources where the accretion rate 
increases and decreases 100 fold during an outburst. Cyg X-1 
stays in the hard state most of the time, but several spectral tran- 
sitions to the soft state were observed. Because of the narrow 
range in mass accretion rate the hard state of Cyg X-1 is still 



intermediate in nature as the observations of Ibragimov et al. 
(2005) show. 

Also by another fact Cyg X-1 differs from the transient 
sources. There is no noticeable hysteresis in Cyg X-1 (Zdziarski 
et al. 2002, Zdziarski & Gierlinski 2005). This was pointed out 
as a problem also by Lachowicz & Czerny (2005). Lightcurves 
of transient sources show a clear hysteresis in the transition lu- 
minosities (the soft/hard transition occurs at luminosities lower 
by a factor of 3 to 5 than the hard/soft transition, e.g. in Aql X-1 
(Maccarone & Coppi 2003), GX 339-4 (Miyamoto et al. 1995, 
Nowak et al. 2002, Zdziarski et al. 2004), XTE J 1650-500 (Rossi 
et al. 2003), XTE J1500-564 (Rodrigues et al. 2003). The differ- 
ent luminosities at spectral transition can be understood as due 
to different maximal evaporation rates in the two states (trigger- 
ing the transition) for either hard or soft irradiation of the corona 
by the central source (Meyer-Hofmeister et al. 2005). 

The new feature discussed in our work here, an inner re- 
condensation disk, offers a solution for the mystery of why Cyg 
X-1 has no hysteresis while all the transients have it. 

In transient sources the re-condensation disk recedes inward 
and finally disappears completely as the mass accretion rate con- 
tinues to fall, as is the case in the decline from outburst into qui- 
escence. Then, no inner disk remains on which re-condensation 
could occur when the system moves into the next outburst. The 
spectrum remains hard at increasing luminosity. Only when the 
accretion rate reaches the higher critical value for which the tran- 
sition under hard irradiation occurs does it become soft. 

Cyg X is different. When that system changes from a soft 
state excursion back to its usual hard state it does so with only 
a modest decrease of the mass accretion rate and never falls into 
deep quiescence. This allows an inner re-condensation disk to 
survive as the observations indicate. When now another episode 
of transition to a soft state occurs the mass accretion rate rises, 
re-condensation increases, the inner disk expands and carries an 
increasing disk mass flow until it smoothly merges with the outer 
disk to form the standard soft accretion state. Together with re- 
condensation the irradiation is completely reversible and no hys- 
teresis occurs. This explanation relates the special character of 
the wind accretor to the otherwise enigmatic feature of "no hys- 
teresis" in a natural way and seems a strong support for the re- 
condensation model. 

10. Discussion 

A more detailed numerical investigation of the re-condensation 
process is desirable. The present simple analytical model might 
serve to show the features and their dependence on the parame- 
ters to be expected in this interesting accretion mode. 

10.1. Interpretation of the intermediate state 

When the mass accretion rate gradually decreases further as a 
source declines from outburst the re-condensation rate onto the 
receding inner disk declines, a significant part of the accretion 
flow stays in the ADAF and shows characteristics of the hard 
state but modified by the existence of cool material in the form 
of a weakly accreting inner disk. From our analysis we find that 
the intermediate states during the outburst decline last longer for 
a slow decrease of the mass flow rate and less time for a fast 
decrease. 

Observations for transient black hole X-ray sources show 
those features. During the slow luminosity decline after the 1999 
outburst of XTE J 1859+226 (Corbel et al. 2004, Fig.9) the spec- 
trum was classified as due to an intermediate state. The steep 
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luminosity decUne of GX 339-4 in 2003 (Remillard 2005, Fig.3) 
is related to the very short intermediate state. 

Short intermediate states of a few days might be caused by 
an inner disk disappearing in its viscous time. Similar to these 
short states is the situation during rise to outburst, where the in- 
ner edge of the disk moves inward within a viscous time, only 
an ADAF inside. If the mass flow rate in the disk is always 
well below the maximal evaporation rate the disk is truncated 
all the time, and no thermal conduction affects the ADAF. XTE 
Jl 1 18H-480 is a well studied source which stays in the hard state 
even during outburst. We would expect a disk truncation at sev- 
eral hundred Rs. But Esin et al. (2001) from multiwavelength 
observations concluded that the disk truncation lies at about 
55/?s in outburst. How can this be understood? The highest mass 
flow rate during outburst was found to be about 0.02 MEdd, the 
value generally accepted for state transition (Maccarone 2003). 
This means that the mass flow rate could have surpassed the 
critical rate for a short while and the outer disk began to ex- 
tend inward. Re-condensation immediately would begin. This 
re-condensation can then continue even as the accretion rate 
drops below the critical rate, leading to an intermediate state. 

10.2. Inner edge of the disk? 

The disk truncation in general varies with the mass accretion 
rate. For low luminosity sources. X-ray binaries or AGN in a 
hard spectral state, the inner edge of the outer disk was deter- 
mined by fitting multiwavelength spectra, understood as caused 
by an ADAF and an outer disk. For very low accretion rates large 
truncation radii were found, of the order of 10^ to lO'^Rs (e.g. 
Esin et al. 1997, Nai-ayan et al. 1997, Narayan 2005). For X-ray 
binaries with higher mass accretion rates in phases of rise to an 
outburst or decay, from a multicolor disk model smaller values 
for the inner disk radius were deduced (compare Sobczak et al. 
2000, see examples in the review of McClintock & Remillard 
2006). In a recent detailed investigation of the spectrum of Cyg 
X-1 Ibragimov et al.(2005) found truncation radii around 50-100 
R5. Note that the mass flow rate was not much below the transi- 
tion rate, i.e. a typical intermediate state. 

Comparing with the predictions of the standard evaporation 
model we find good agreement for the low luminosity sources. 
The model also is consistent with no disk truncation for high 
mass flow rates. 

For the intermediate state of Cyg X-1, Ibragimov et al. 
(2005) analyzed spectral observations in terms of Compton re- 
flection, the index of the hard power law spectrum, the width 
of the Fe Ka line, and an additional soft excess of a few keV 
of unknown origin. With the correlation between photon index 
and amplitude of reflection they inferred disk cut-off' distances 
of 50-100 R5, varying with the accretion rate. Esin et al. (2001) 
inferred a similar truncation at 55 Rs for XTE Jl 1 18-1-480 in an 
apparently similar hard (-intermediate) state (see the discussion 
in the preceding subsection). 

On the other hand. Miller et al. (2006a) argue for an un- 
truncated disk extending to the last stable orbit during the low- 
hard state of the 2004 outburst of GX 339-4, at a luminosity of 
5% of the Eddington value. From a preliminary analysis of data 
from SWIFT J17335-0127 in the low-hard state of its 2005 out- 
burst. Miller et al. (2006b) see indication for a similarly untrun- 
cated disk, here at a luminosity of one third of a per cent of the 
Eddington value. (Note that in optical thin emission the accre- 
tion luminosity can be lower than the corresponding accretion 
rate, both measured in their Eddington units.) This would sug- 



gest an untruncated re-condensation disk in these systems. The 
situation deserves further clarification. 

The re-condensation model discussed here leaves open the 
question of if and where a remaining inner disk may itself be 
truncated. Neglected terms in the analysis, different assumptions 
about y, and a more detailed calculation all can affect the out- 
come of the re-condensation rate estimates. One also cannot ex- 
clude that a parameter used to describe the ADAF varies with 
distance. The ratio of magnetic pressure to gas pressure, for ex- 
ample, could increase at smaller radii as accretion carries mag- 
netic flux inwards. This decreases /3 and C (Eq.dTTTl). If C falls 
below one (Eq.(l23]l). condensation turns into evaporation and 
the disk can become truncated before the final last stable orbit 
is reached. The effect of disk magnetic fields on evaporation in 
the standard model is discussed in Meyer & Meyer-Hofmeister 
(2002). 

10.3. Intermediate spectral states in disks around 
supermassive black holes? 

The solutions discussed in this paper are independent of the mass 
m of the black hole. Thus they should be applicable also to ac- 
cretion on supermassive black holes in AGN. Recently Jester 
(2005) studied the distribution of AGN bolometric luminosi- 
ties and black hole masses for objects from the SDSS spectro- 
scopic quasar survey to test the existence of two different accre- 
tion modes (as predicted theoretically, Narayan et al. 1998) and 
found a change of mode at an Eddington-scaled accretion rate of 
about 0.01 . Markowitz and Uttley (2005), comparing power den- 
sity functions pointed out the analogy of low luminosity AGN 
to low/hard state black hole X-ray binaries. This suggests that 
some of the AGN accretors should be found in an intermedi- 
ate state like their galactic stellar mass counterparts. To actually 
observe a transition through the different accretion states might 
be difficult as dynamical and viscous timescales scale with the 
mass of the black hole. Yuan et al. (2004) compared ROSAT and 
XMM Newton observations of 386 sources and found in one of 
them, the Seyfert-LINER galaxy NGC 7589, an increase in X- 
ray flux by a factor of >10 over a time of 5 years, a timescale of 
the order of the diffusion time at a distance of a few hundred 
Schwarzschild radii from a 10^ Mq black hole. As the XMM 
high state luminosity was estimated to be a few percent of the 
Eddington value the authors suggest that the system might have 
been found during a spectral transition. 

11. Conclusions 

We have presented a simple analytical model for re-condensation 
of gas from an ADAF onto a disk below under the action of 
thermal conduction and radiation. The model is able to explain 
how cool gas can exist in the close neighborhood of an accret- 
ing black hole, after a standard accretion disk has become trun- 
cated and an ADAF has formed, a situation that occurs when the 
mass flow rate decreases during the transition from soft to hard 
spectral state of galactic black hole X-ray binaries. A fairly com- 
plete re-condensation of the ADAF into the disk would extend 
an inner thin disk accretion with its soft spectral state to lower 
mass accretion rates, and suggests a hysteresis effect indepen- 
dent of and in addition to the irradiation effect discussed ear- 
lier. The model also clarifies why in soft state the considerable 
coronal flow at distances of several hundred Schwarzschild radii 
predicted theoretically does not continue as a hot 2-temperature 
flow to the interior. Such a flow would appear as a hard contri- 
bution in the spectrum that is not observed. The model receives 
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considerable support from being able to resolve the mystery of 
why Cyg X-1 does not show hysteresis between the luminosity 
at soft/hard and hard/soft spectral transition that all X-ray tran- 
sient sources display. 

This simple model deserves further, more detailed investiga- 
tion. 
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